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ABSTRACT 



Context. Masses are a fundamental parameter, but they are not well known for most hot subdwarfs. In general, the mass of a hot 
subdwarf is derived with asteroseismology or dynamical methods, for which it is often difficult to obtain the necessary data from 
observations. 

Aims. We intend to find an approach to deriving the masses of hot subdwarfs from observational data in the literature. 

Methods. We presented full evolutionary calculations for hot subdwarfs in a wide mass range (0.33 Mq to 1.4 Mq) for a Population I 

J.4 

metallicity of Z=0.02, and obtained a relation between Mp and log(-^), where Mp, T^-g, and g axe. the most probable mass, effective 
temperature, and gravity. This relation is used to study the masses of some observed hot subdwarfs. 

Results. We proposed a method of determining the masses of hot subdwarfs. Using this method, we studied the masses of hot 
subdwarfs from the ESO supernova la progenitor survey and Hamburg quasar survey. The study shows that most of subdwarf B stars 
have masses between 0.42 and 0.54 Mq, whilst most sdO stars are in the range 0.40 ~ 0.55 Mq. Comparing our study to the theoretical 
mass distributions of Han et al. (2003), we found that sdO stars with mass less than ~ 0.5 Mq may evolve from sdB stars, whilst most 
high-mass(> 0.5 Mq) sdO stars result from mergers directly. 

Key words, stars: fundamental parameters - stars: horizontal-branch - subdwarfs 
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1. Introduction 

Hot subdwarfs are generally classfied into three types by their 
spectra. These are subdwarf B (sdB, with a surface effective tem- 
perature, T^ff, in a range from 20,000 to 40,000K, with H-Balmer 
absorption lines wider than in normal B stars), subdwarf O (sdO, 
Teff ranging from 40,000K to 80,000K with strong He absorp- 
tion lines), and subdwarf OB (sdOB, a transition between O and 
B, Moehler et al. 1990; Heber 2009). These objects are located 
below the upper main sequence on the Hertzsprung-Russell di- 
agram(HRD) and are also known as extreme horizontal branch 
(EHB) stars from the view of their evolutionary stages; i.e., they 
are believed to be core He-burning objects with extremely thin 
hydrogen envelopes (< 0.02 Mq). 

Hot subdwarfs are an important population in several re- 
spects. For example, pulsating sdB stars are standard candles in 
distance determination (Kilkenny et al. 1999). Likewise, close 
binaries composed of an sdB star and a massive white dwarf 
(WD) are qualified as supernova la progenitors (Maxted et al. 

2000) . Moreover, hot subdwarfs are an important source of far- 
UV light in the galaxy, and they are successfully used to explain 
the UV-upturn in elliptical galaxies (Kilkenny et al. 1997; Han 
et al. 2007). 

Several scenarios have been proposed to explain the forma- 
tion of these objects, i.e. strong mass loss for a star on the red 
giant branch (Dorman et al. 1993, D'Cruz et al. 1996), mass 
transfer in close binary systems (Mengel et al. 1976), and the 
coalescence of two helium white dwarfs (He-WDs) (Iben 1990, 
Webbink 1984). Recent radial velocity (RV) surveys reveal that 
two thirds of all sdB stars reside in close binaries (Maxted et al. 

2001) . Han et al. (2002, 2003) propose a binary model for the 



formation of hot subdwarf stars, in which three channels (sta- 
ble Roche lobe overflow (RLOF), common envelope (CE) ejec- 
tion, and merging of two He-WDs) are included. (For details see 
Han et al. 2002, 2003.) This binary model successfully explained 
most of the observational characteristics of hot subdwarfs, and is 
now widely used in the study of hot subdwarf stars, e.g. O'Toole, 
Heber & Benjamin, 2004. 

Mass is a fundamental parameter of stars, but it is very un- 
certain for most hot subdwarfs, which are assumed to be around 
0.5 Mq (Heber 1986, Saffer et al. 1994). Sometimes we need a 
more precise mass to study these objects further In recent years, 
more than 2300 objects have been included in a hand hot sub- 
dwarf database (0stensen 2004), among these more than 200 
hot subdwarfs have been studied in detail for atmospheric pa- 
rameters (e.g. Saffer et al. 1994; Maxted et al. 2001; Edelmann 
et al. 2003, Lisker et al. 2005, Stroeer et al. 2007). However, 
only a few hot subdwarfs have well-defined masses from as- 
teroseismology, e.g. PG 0014+067 (Brassard et al. 2001), PG 
1219-H534 (Charpinet et al. 2005a), PG 1325-Hl01(Charpinet et 
al. 2006), and Feige 48 (Charpinet et al. 2005b). In addition. 
Wood & Saffer (1999) show an sdB star HW Vir (PG 1241-082), 
which is in a double-lined eclipsing binary, to be 0.48 ± O.O9M0 
from dynamical methods. For other hot subdwarfs, which are not 
pulsating stars or are in binary systems without detailed light or 
radial-velocity curves, masses are derived from theoretical evo- 
lutionary tracks (Dorman et al. 1993), e.g., choosing the track 
closest to the data point on the T^ff - log(^) diagram, where Teff 
and g are effective temperature and surface gravity, respectively. 
The commonly used evolution library of hot subdwarfs is that of 
Dorman et al. (1993), with masses from 0.4 to 0.5 Mq. However, 
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Fig. 1. The maximum envelope masses for different core masses 
in our hot subdwarf models (squares). The solid line is a linear 
fitting of the square points, with errors less than 0.00 IMq. 

Han et al.(2002, 2003) show that the masses of hot subdwarfs 
can be in a much wider range; i.e., from 0.3 Mq to more than 0.7 
Mq. Thus, it is necessary to extend the evolution library of hot 
subdwarfs or find a simple method of estimating the masses of 
hot subdwarfs. 

In this letter, we carry out full evolutionary calculations for 
hot subdwarfs and obtain an approach to determining the masses 
of hot subdwarfs. Using this approach, we study the masses and 
mass distributions of observed hot subdwarfs. 



2. The model 

The stellar evolution code used was originally developed by 
Eggleton (1971, 1972, 1973). The code has been updated with 
the latest input physics over the past three decades as described 
by Han, Podsiadlowski & Eggleton (1994, hereafter HPE) and 
by Pols et al. (1995, 1998). We set the ratio of mixing length to 
local pressure scale height, a = l/Hp, to 2.0 and the convective 
overshooting parameter, (Jqv, to 0.12 (Pols et al. 1997; Schroder 
et al. 1997). The opacity tables for various metallicities are com- 
piled by Chen & Tout (2007) from Iglesias & Rogers (1996) and 
Alexander & Ferguson (1994). The initial hydrogen mass frac- 
tion, X, is obtained by, X - 0.76-3.0Z, where Z is the metallicity 
(Pols et al. 1998). 

We constructed a series of zero-age extreme horizontal 
branch models (ZAEHB, i.e. helium begins stably burning in the 
core) of Pop I (Z=0.02) and carried out their evolutions in detail. 
The model grid contains a wide mass range; i.e., the core mass 
ranges from 0.33 to 1.4 Mq, and the envelope mass changes, by 
a step of about 0.005 Mq, from zero to the maximum envelope 
mass (see Fig[T]), beyond which the constructed object is not a 
hot subdwarf any more (i.e. T^ff < 20,000K or log(^) < 4.5). 

As an example, we show the evolutionary tracks of hot sub- 
dwarfs with a core mass of 0.475 Mq in Fig. 2. The envelope 
masses, Menv, are 0.000, 0.001, 0.005, 0.010, and 0.015 Mq. 
TAEHB is the termination-age EHB, i.e., at the core helium ex- 
haustion (Dorman et al. 1993). The region between ZAEHB and 
TAEHB is known as the main sequence. The age difference be- 
tween two neighbouring crosses is 10 Myr. From the figure, we 
see that most of the observed sdB stars are located in the main se- 
quence. However, we cannot confirm their core masses and their 
evolutionary stages yet, since a star with a different core mass in 
a different evolutionary stage is also likely in this region of the 
Teff - log(§) diagram (see Sect.3). For instance, an sdO star is 
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Fig. 2. Evolutionary tracks of hot subdwarfs with a core of 0.475 
Mq and Z=0.02. The envelope masses are 0.000, 0.001, 0.005, 
0.010, and 0.015 Mq. ZAEHB and TAEHB represent the on- 
set and the end of main sequence, respectively (see the text for 
detail). The circles and diamonds are for sdB and sdO stars, re- 
spectively, from the ESO supernova la progenitor survey (SPY, 
Lisker et al. 2005, Stroeer et al. 2007). The age difference be- 
tween two neighbouring crosses is 10 Myr 

possibly a main sequence star with a core mass higher than that 
of an sdB star, or in the stage of post main sequence of a sdB 
star. 

3. Methods 

As we know, more than one evolutionary tracks pass through 
the same point on the T^ff - log(g) diagram. Thus, we cannot 
figure out the most probable one from them just by comparing 
the evolutionary tracks. To solve this problem, we intend to find 
another approach to deriving the masses of hot subdwarfs. 

Hot subdwarfs are core helium-burning stars with hydro- 
gen envelopes that are too thin to sustain hydrogen burning; 
therefore, the luminosity of a hot subdwarf, L, depends upon 
the core mass, M^, and the age, t. Meanwhile, we know that 

L ~ T\ ■ R- ~ Mr ■ — , where R is the radius. Thus, we have 

eff * 4 

a relation as f{M^, t) ~ log(^), where Teff and g are two fun- 
damental parameters of hot subdwarfs from observations. Thus, 
we choose log(^) as a basic parameter for further study of the 
masses of these objects. 

In general, stars evolves at different speeds in different evo- 
lutionary stages. If a star evolves slowly in a certain stage, it 
has more chance of being observed in this stage. We may then 
find out the most probable mass via studying the change rate 

of log(-^), as described in the following. We divide log(-ji) 
into small intervals of 0.1 from 10.0 to 16.0, and calculate the 

evolutionary time, 6t, of each model in these intervals. Thus, for 

7.4 

each interval of log(-|ii), we may obtain a core mass, with which 
a hot subdwarf stays in it for the longest time amongst all the 
models (i.e. having the maximum value of 6t), and has the most 
chance of being observed. This core mass is then the most prob- 
able mass, Mp, in the observations, and is assumed to be the real 
mass of the subdwarfQ We finally obtained a relation between 

y4 

Mp and log(-^), which is shown in Fig. 3. From this relation we 

' Note that the mass given here is the core mass, which is a lower 
limit of the total mass of hot subdwarfs. 
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Fig. 3. The Mp-log(^) relation. Mp is the most probable mass. 



may simply know the mass of a hot subdwarf for given Teg and 
8- 

Using the Mp-log(-^) relation, we studied the masses of 
some observed hot subdwarfs, including 164 sdB stars and 46 
sdO stars. Some sdB stars are from the ESO supernova la pro- 
genitor survey (SPY, Lisker et al. 2005) and the other sdB 
stars are from the Hamburg quasar survey (HQS, Edelmann et 
al.2003). The sdO stars come from SPY (Stroeer et al. 2007). 
The mass distributions for sdB and sdO stars are shown in Figs .4 
and 5, respectively. 

The error on the mass, AMp, depends on the parameter 

log(^), and this dependence is shown in the upper panel of 

Fig. 6, where f (= -rr—^hi—J is derived from the MD-log(— ) 

relation. The log(-^) distributions of sdB and sdO stars from 
the SPY and the HQS are shown in the bottom panel of this fig- 
ure. From this figure we see that / is lower than 0.5 for most 
hot subdwarfs. To discuss the error on the mass obtained from 
our method, we should consider the error from both theory and 

observation. The Mp-log{^) relation here is derived from Pop 
I models (Z=0.02). Thus, the different metalicities result in a 
different log(^) for the same mass. For example, Alog(^) 
is about 0.04 dex for the models between Z - 0.0001 and Z 
= 0.06. For stars from the ESO supernova la progenitor sur- 
vey, the observational errors for sdB stars are ATeff - 36QK and 
A log(g) - 0.05 dex, and the observational errors for sdO stars 
are Alog(reff) - 0.011 dex and Alog(g) = 0.097 dex, respec- 
tively. These correspond to Alog{T*g/g) of 0.07 dex and 0.14 
dex. Thus, the total errors, Alog(T^g/g) (theory plus observa- 
tion), are 0.11 dex and 0.18 dex, respectively. According to the 

/ ~ log(^) relation shown in Fig. 6, the errors on the masses 
are less than 0.055 Mq and 0.09 Mq for most sdB and sdO stars, 
respectively. For sdB stars from the Hamburg Quasar Survey, we 
derived the mass errors for each star in a similar way. The results 
show that the majority of errors are in the 0.097 M© to 0. 1 17 Mq 



Fig. 4. The mass distribution of sdB stars (derived from the Mp- 

log(-^) relation) from the ESO supernova la progenitor survey 
(SPY, Lisker et al. 2005) and the Hamburg quasar survey (HQS, 
Edelmann et al.2003). As a comparison, we also show the mass 
distribution from the theoretical study of Han et al. (2003) and 
from asteroseismology (Fontaine et al. 2006). 

range(95 percent confidence interval), while the mean error is 
0.107 Mq. 

4. Results and discussion 

As shown in Fig. 4, the sdB stars from SPY and HQS have the 
same mass distributions — the KS test shows that the two samples 
come from the same distribution at a high level of confidence(> 
93.5%). Most of the sdB stars have masses ranging from 0.42 
to 0.54 Mq, and the mean mass is about 0.50 Mq, equal to that 
assumed in previous studies. As a comparison, we also show the 
mass distributions from the theoretical study of Han et al. (2003) 
0and from asteroseismology (Fontaine et al. 2006) in Fig. 4, from 
which we see that the distribution obtained in this paper matches 
the other two closely. 

Figure 5 shows that most sdO stars are in 0.40 ~ 0.55 Mq, 
although some sdO stars have much higher masses than these 
values. Since low-mass hot subdwarfs (e.g. less than ~0.5 Mq) 
can not reach 40,000K (the minimum temperature for sdO stars) 
during the main sequence (see Fig. 2 in Han et al. 2002), sdO 
stars with mass less than ~0.5 Mq most likely evolve from sdB 
stars, i.e., these objects initially appear as sdB stars when they 
are on the main sequence, and as sdO stars when they evolved 
off the main sequence (post-EHB). The typical evolutionary time 
scales for hot subdwarfs on the MS and post-EHB are ~ 160 and 
~ 20Myr, respectively. Thus, observationally the number of sdO 
stars (evolving from sdB stars) is expected to be about 12% of 
that of sdB stars, and for SPY this equates to about 9 objects. 
They are likely He-deficient sdO stars. There are 13 He-deficient 
sdO stars in SPY, and 7 ones (likely evolve from sdB stars) of 

- In the paper of Han et al. (2003), the mass distribution are for both 
sdB and sdO stars. We seperated these objects, thus the distribution pre- 
sented here is different from that of Han et al.(2003). 



4 



X. Zhang et al.: The masses of hot subdwarfs 





SdB-SPY 


H 


sdB-HQS 


□ 


sdO-SPY 




Fig. 5. The mass distribution of sdO stars (derived from the Mp- 

log(-^) relation) from the ESO supernova la progenitor survey 
(SPY, Lisker et al. 2005). The mass distributions from the CE 
and merger scenarios (Han et al. 2003)^ are also shown in this 
figure. 



them with masses less than ~0.5 in our study, similar to the 
9 predicted by theory. 

To more easily compare with theoretical results, we pre- 
sented the mass distributions of sdO stars obtained from Han 
et al. (2003)^ in Fig.5, where the dot-dashed line is from CE 
ejection and the circle-line is from the merger of two He-WDs. 
We see that CE ejection can only account for sdO stars with 
masses lower than ~0.5 Mq. The He-He merger mainly con- 
tributes to high-mass sdO stars, but it is also responsible for a 
small fraction of low-mass sdO stars. In particular, the fraction 
of very low-mass stars (< 0.4Mq) is smaller among sdO stars 
than among sdB stars, as shown in Figs. 4 and 5. This difference 
probably stems from very few sdB stars with masses lower than 

0. 4M0 being hotter than 40,000K during their post-EHB stage. 

We know that a binary resulting from CE ejection gener- 
ally has a short orbital period. Thus, if low-mass sdO stars are 
mainly from the CE ejection channel, then a large fraction of 
these objects should be observed in short orbital period bina- 
ries. Among 23 sdO stars with masses less than ~0.5 M©, (if 
they indeed evolve from sdB stars), we would expect a fraction 
of binaries similar to that of sdB stars, i.e. 39% (Napiwotzki 
et al. 2004). However, only five sdO stars have been identified 
in binaries from SPY. The possible reasons for the low fraction 
are (a) the He-He merger contributes a small fraction of low- 
mass sdO stars, which are all singles; (b)low-mass sdO stars 
from CE ejection generally have lower mass companions (most 
likely M, see Fig. 15 in Han et al. 2003), and may be identified 
as singles; (c) some other observational effects, such as random 
inclination of orbital planes, contribute to a decrease in detec- 
tion efficiency. Moreover, in comparison to the merger channel, 
CE ejection, which may leave a thin hydrogen-rich envelope af- 
ter the ejection, often produces hot subdwarfs with surface he- 
lium deficiency. Thus, if sdO stars are identified in binaries with 
short orbital periods, they are most likely to be helium deficient, 
since the merger of He-He WDs can only produce single helium- 
enriched sdO stars. The observations confirm this prediction; 

1. e., four of the five identified binary sdO stars have short or- 
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Fig. 6. The upper panel: the dependence of the mass errors 
on log(^), derived from the Mp-log(^) relation, where / = 

AM 

., ,Ji I The bottom panel: the log(— ) distributions of sdB 
and sdO stars from the SPY and the HQS. 



bital periods and are helium deficient. An extreme case from the 
CE ejection channel is that all of the common envelope has been 
ejected, leading to a naked He core (hot subdwarfs) and a com- 
panion. This may produce helium-enriched sdO stars with short 
orbital periods. Since this case is rare, the fraction of helium en- 
riched sdO binaries among He-sdO stars (most of them are from 
the merger channel) is very low (4% at most, Napiwotzki et al. 
2004, see also Heber 2008). 
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ABSTRACT 



Context. Masses are a fundamental parameter, but they are not well known for most hot subdwarfs. In general, the mass of a hot 
subdwarf is derived with asteroseismology or dynamical methods, for which it is often difficult to obtain the necessary data from 
observations. 

Aims. We intend to find an approach to deriving the masses of hot subdwarfs from observational data in the literature. 

Methods. We presented full evolutionary calculations for hot subdwarfs in a wide mass range (0.33 Mq to 1.4 Mq) for a Population I 

J.4 

metallicity of Z=0.02, and obtained a relation between Mp and log(-^), where Mp, T^-g, and g axe. the most probable mass, effective 
temperature, and gravity. This relation is used to study the masses of some observed hot subdwarfs. 

Results. We proposed a method of determining the masses of hot subdwarfs. Using this method, we studied the masses of hot 
subdwarfs from the ESO supernova la progenitor survey and Hamburg quasar survey. The study shows that most of subdwarf B stars 
have masses between 0.42 and 0.54 Mq, whilst most sdO stars are in the range 0.40 ~ 0.55 Mq. Comparing our study to the theoretical 
mass distributions of Han et al. (2003), we found that sdO stars with mass less than ~ 0.5 Mq may evolve from sdB stars, whilst most 
high-mass(> 0.5 Mq) sdO stars result from mergers directly. 
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1. Introduction 

Hot subdwarfs are generally classfied into three types by their 
spectra. These are subdwarf B (sdB, with a surface effective tem- 
perature, T^ff, in a range from 20,000 to 40,000K, with H-Balmer 
absorption lines wider than in normal B stars), subdwarf O (sdO, 
Teff ranging from 40,000K to 80,000K with strong He absorp- 
tion lines), and subdwarf OB (sdOB, a transition between O and 
B, Moehler et al. 1990; Heber 2009). These objects are located 
below the upper main sequence on the Hertzsprung-Russell di- 
agram(HRD) and are also known as extreme horizontal branch 
(EHB) stars from the view of their evolutionary stages; i.e., they 
are believed to be core He-burning objects with extremely thin 
hydrogen envelopes (< 0.02 Mq). 

Hot subdwarfs are an important population in several re- 
spects. For example, pulsating sdB stars are standard candles in 
distance determination (Kilkenny et al. 1999). Likewise, close 
binaries composed of an sdB star and a massive white dwarf 
(WD) are qualified as supernova la progenitors (Maxted et al. 

2000) . Moreover, hot subdwarfs are an important source of far- 
UV light in the galaxy, and they are successfully used to explain 
the UV-upturn in elliptical galaxies (Kilkenny et al. 1997; Han 
et al. 2007). 

Several scenarios have been proposed to explain the forma- 
tion of these objects, i.e. strong mass loss for a star on the red 
giant branch (Dorman et al. 1993, D'Cruz et al. 1996), mass 
transfer in close binary systems (Mengel et al. 1976), and the 
coalescence of two helium white dwarfs (He-WDs) (Iben 1990, 
Webbink 1984). Recent radial velocity (RV) surveys reveal that 
two thirds of all sdB stars reside in close binaries (Maxted et al. 

2001) . Han et al. (2002, 2003) propose a binary model for the 



formation of hot subdwarf stars, in which three channels (sta- 
ble Roche lobe overflow (RLOF), common envelope (CE) ejec- 
tion, and merging of two He-WDs) are included. (For details see 
Han et al. 2002, 2003.) This binary model successfully explained 
most of the observational characteristics of hot subdwarfs, and is 
now widely used in the study of hot subdwarf stars, e.g. O'Toole, 
Heber & Benjamin, 2004. 

Mass is a fundamental parameter of stars, but it is very un- 
certain for most hot subdwarfs, which are assumed to be around 
0.5 Mq (Heber 1986, Saffer et al. 1994). Sometimes we need a 
more precise mass to study these objects further In recent years, 
more than 2300 objects have been included in a hand hot sub- 
dwarf database (0stensen 2004), among these more than 200 
hot subdwarfs have been studied in detail for atmospheric pa- 
rameters (e.g. Saffer et al. 1994; Maxted et al. 2001; Edelmann 
et al. 2003, Lisker et al. 2005, Stroeer et al. 2007). However, 
only a few hot subdwarfs have well-defined masses from as- 
teroseismology, e.g. PG 0014+067 (Brassard et al. 2001), PG 
1219-H534 (Charpinet et al. 2005a), PG 1325-Hl01(Charpinet et 
al. 2006), and Feige 48 (Charpinet et al. 2005b). In addition. 
Wood & Saffer (1999) show an sdB star HW Vir (PG 1241-082), 
which is in a double-lined eclipsing binary, to be 0.48 ± O.O9M0 
from dynamical methods. For other hot subdwarfs, which are not 
pulsating stars or are in binary systems without detailed light or 
radial-velocity curves, masses are derived from theoretical evo- 
lutionary tracks (Dorman et al. 1993), e.g., choosing the track 
closest to the data point on the T^ff - log(^) diagram, where Teff 
and g are effective temperature and surface gravity, respectively. 
The commonly used evolution library of hot subdwarfs is that of 
Dorman et al. (1993), with masses from 0.4 to 0.5 Mq. However, 



